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Abstract 

The non isotropic and non central elliptical shape distributions via the Le and 
Kendall SVD decomposition approach are derived in this paper in the context of 
invariant polynomials and zonal polynomials. The so termed cone and disk densities 
here obtained generalise some results of the literature. Finally, some particular 
densities are applied in a classical data of Biology, and the inference is performed 
after choosing the best model by using a modified BIC criterion. 

1 Introduction and the main principle. 

The multivariate statistical shape theory has been develo ped in the last two decades ar ound 
the classical works based o n normality and isotropy, see iGoodall and Mardia I (Il993h and 



Drvden and Mardia I (ll998h and the references there in. Recent works extended this results 



to elliptical model s and partial non isotropy, see for example Caro-Lopera et al ( 20091) and 



iDiaz-Garcia et al. ( 2003 ). however some important problems remain, the study of the shape 
theory without any restriction of the covariance matrix in the elliptical model. 

The problem arises from the point of view of applications, the isotropic assumption 
= Ik for an elliptical shape model of the form 

X ~£Ar X A'(At x ,£ x ,0,/j), 

restricts substantially the correlations of the landmarks in the figure. Then, we expect the 
non isotropic model, with any positive definite matrix 0, as the best model for considering 
all the possible correlations among the anatomical (geometrical o mathematical) points. 
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This problem can be solved by considering the following procedure: Let be 



X ~ £nxk(h x , S x , ©, h), 
if 1 / 2 is the positive definite sq uare root of the matrix 0, i .e. = (0 1 / 2 ) 2 , with 1 / 2 : 



K x K . iGupta and Vargal (|1993l p. 11), and noting that 



where 
then 



X0 ] X' = X(0" 1/2 0- 1/2 )- 1 X' = X0- 1/2 (X0~ 1/2 )' = ZZ', 

Z = X©~ 1/2 , 



Z ~ fAr X Ar(At z ,5] x ,I^,/l) 

with n z = /i x 0" 1/2 , fsee lGupta and Vargal (|l993l p. 20)). 



And we arrive at the classical starting point in shape theory where the original landmark 
matrix is replaced by Z = Then we can proceed as usual, removing from Z, 

translation, scale, rotation and/or reflection in order to obtain the shape of Z (or X) via the 
QR, SVD, polar decompositions. The shape theory as sociated with t he SVD decomposition 
ca n be study from two diffe rent approaches, one due to lGoodalll (|199ll ) and another proposed 
bv lLe and Kendall I (Il993l) . 



We study in this paper the statistical approach of Le an d Kendall under a gener alised 
elliptical model. First, recall some facts of this technique (|Le and Kendall ( 1993h ). It 



is known that the shape of an object is all geometrical information that remains after 
filtering out translation, rotation and scale information of an original figure (represented 
by a matrix X) comprised in N landmarks in K dimensions. Hence, we say that two 
figures, Xi : N X K and X2 : N x K have the same shape if they are related with a 
special similarity transformation X 2 = (3X. 1 H + InI', where H : K x K G SO(K) (the 
rotation), 7 : Jf x 1 (the translation), ljy : N x 1, ljy = (1, 1, . . . , 1)', and f3 > (the scale). 
Thus, in this context, the shape of a matrix X is all the geometrical information about X 
that is invariant under Euclidean similarity transformations. Then, the shape space is the 
set of all possible shapes, it is the orbit space of the non-coincident N landmarks in $l K 
under the action of the Euclidean similarity transformations. The dimension of this space is 
NK — K — 1 — K(K — l)/2, it is, the original dimension NK is reduced by K for location, 
by 1 for uniform scale and by K(K — l)/2 for rotation. In other words, the shape of X is the 
set {PT : T e SO{K)} where P is the so termed pre-shape of X defined as P = LX/||LX|| 
(L is Helmert submatrix, for example) which is invariant under translation and scaling of X. 
The rotated P on the pre-shape sphere is termed a fibre of the pre-shape space S^, these 
fibres do not overlap and corresponds one to one with shapes in the shape space E^, it is, 
the pre-shape space is partitioned into fibres by the rotation group SO(K) and the fibre is 
the orbit of P under the action of SO(K). Thus, is the quotient space of under the 
action of SO(K), in notation = / SO(K), which means that the shape of X is an 
equivalent class under the action of the group of similarity transformations. 

Now, the statistical theory of shape associated to this approach studies the effect of 
randomness and assume a probabilistic model for the original matrix in order to obtain the 
density of the pre-shape (cone) and shape (disk). The complete procedure for obtaining the 
shape of an original X can be summarised in the following steps: 

LX0 1/2 = LZ = Y = V'DH = rV'WH = ? V'W(u)H, 

where the matrix L : (N — 1) x N has orthonormal rows to 1 = (1, . . . , 1)'. L can be a 
submatrix of the Helmert matrix, for example. Here Y = V'DH is the SVD of matrix 



2 



Y, with V 
diag(L>!, . . . 



n X (N - 1) and H 
D n ); W = D/r, r = | 



D 



n x K semiorthogonal matrices and D 



n x n. D 



Y 



Then the standard problem considers a model for X and finds the so termed cone and 
disk densities, which are the densities of D and W(u), respectively. A number of published 
works have studied the classical problem which assumes a non central Gaussian model for 
X very restricted by the isotropy assumption = 1^- and Ex = <j 2 In in 

X~7V JV x^(M x ,S x ,0). 

These restrictions facilitates the integration in terms of the zonal polynomials and the asymp- 
totic densities can be derived and applied. This procedure based on normality and isotropy 
is very common in literatur e of shape under QR, SVD and affine decompositions too, see 
Goodall and Mardia I ([19931 ) and iDrvden and Mardia I (|1998l ). and the references there in. 

However, it is clear that the gaussian case do not support all the applications and the 
statistical theory of shape could enriched if complete families of cone and disk densities 
are available for a particular experiment and the re searcher c a n mo d el the situation by 
applyi n g a model sele ctio n criteria (see for exam ple iRissanenl (|l978l ). H ass and Rafterv 
(ll995h . [Rafte^v1 (Il995h and lYang and Yand (l2007h . among many others). 



Therefore, in this paper we propose the statistical theory of shape of Le and Kendall's 
approach under any non central elliptical model without any restriction of the covariance 
matrix. 

Explicitly, consider a full covariance elliptical model, indexed by the generator function 

h, 

X ~ £Ar x K(/X x ,S x ,0,/l), 

so, by the main principle given above, at the beginning of this section, we have that 

Z ~ £jvxif (/^z, S x , 1 K , h), 

with Z = X0~ 1/2 , /x z = // x © 1/2 . 

Then the (SVD) Le and Kendall's shape coordinates u of X arc constructed in several 
steps summarised in the expression 

LX0 1/2 = LZ = Y = V'DH = rV'WH = rV'W(u)H. (1) 

Denote = L/x x , so Y : (N — 1) x K is invariant to translations of the figure Z, and 

Y ~ ^_ lxK (^0- 1/2 , S ® I K , h), 

where X = LSxL', meanwhile the matrix W, the shape of X, is invariant under (transla- 
tions), rotations and scaling of the landmark data matrix X. 

By considering the above main principe, this paper studies Le and Kendall's approach for 
the shape theory based on the SVD decomposition and any non isotropic non central elliptical 
model. Section [5] obtains the general densities, the so termed cone and disk densities, with 
some corollaries. Then, the central case and its invariance is studied in Section[3] At the end 
inference on small and large mouse vertebra data is performed with the classical Gaussian 
model and two non normal Kotz models, then the best model is chosen by a modified BIC 
criterion and the corresponding test for equality in mean disk shape is obtained. 



2 Shape Theory via SVD Le and Kendall's approach 

We start with the jacobian of the corresponding decomposition. 
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Lemma 2.1. Let be Y : N — 1 x K , then there exist V 6 V n _N-i represents the Stiefel 
manifold, H G V n ,K and D : n x n, D = diag(Di, . . . ,D n ), n = min[(iV — l),K]; D\ > 
D2 > ■ ■ ■ > D n > 0, such that Y = V'DH; This factorisation is termed non-singular part 
of the SVD. Then 



(dY) = 2~™|D 



[](A 2 - Dj)(dT>)(VdV')(HdH'). 



Proof. See lDiaz-Garcia et all (|l997f) . □ 
In order to o btain the joint density function of (V, D) we need the following generalisation 
of l.lamesl (jl964l eq. (22)). 

Lemma 2.2. Let X : K X n, Y : K x K and H 6 V n ,K- Then 
1. 

[ [tr(Y + XUmUdW) = ^ V V ^ ^^ XX ') 



(^)a 



where |(tr Y)^ 1 tr XH < 1 and trY^O. 
/ tr(Y + XH)etr(r-(Y + XH))(HrfH') 

JH£V„,K 



2 n JT Kn/2 



1 „ 



T n [\K] 



— etr(rY) tr Y F l( -^; -XX') + E E TI 



(/ + ±)C A (±XX) 



/=o a ^2 



where p £ I, r £ $, C K (B) are i/ie zonal polynomials of B corresponding to the partition 

« = (a, •••/?) 0/ /, w^/i x)i=i /« = /' arirf ( a )« = n l= i( a - u - i )/ 2 )f j > ( a )f = a ( a + 

1) ■ ■ ■ ( a +f— l), are the generalised hypergeometric coefficients and^Fi is the Bess el function, 
Jamei \l96i ). 



Proof. 



1. From Lemma 9.5.3 lMuirheadl (|1982l Lemma 9.5.3, p. 397) we have 

2 n Kn/I 

lH-Y„ ■■; F n [^K} JO(K) 



r 2™ TT Kn / 2 r 

/ [tr(Y + xh)F (mm') = — _ / 



[tr(Y + XH)] p (dH). 



Furthermore, for tr Y ^ and |(tr Y) 1 trXH| < 1 



[tr(Y + XH)f = (trY) P E^T"( trY ) _/ ( trXH ) 



/! 



Now from I James! (|l964 eqs. (46) and (22))) it follows that 
/ H6V [tr{Y + XH)]*(HdH') = 



2 n 7T Kn/2 (p) 2/ (trY)- 2 / (j)/ Ca(xx0- 



A 



(2/)! {\K) X 



the result follows, noting that (i) / /(2/)! = 1/(4^/!) and that C A (aXX') = a- f C* A (XX'). 
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2. This follows b y expanding the exponentials in series of powers and by applying (22) 
and (27) from 1 James! (|l964h . 

□ 

Thus, we can obtain: 
Theorem 2.1. The joint density of (V, D) is 

n nf |D|7V _ 1+ K- 2ll JJ(£>2 _ D 2) 

/v,d(V,D) = 



x E E 1 \ U i K) (^5T 1 V'D 2 V) . 

t=o k t -V2"; K 

Proof. Let be 17 = X _1 /i.0 _1 /i', then the density of Y is given by 
/ Y (Y) = —^h [tr (S- 1 YY / + Si) - 2tr/x'2T x Y] . 

|S| 2 

Now, make the change of variables Y = V'DH, so, by Lemma I2TT1 the joint density function 
of V, D, H is 



2 -n| D |JV-l+*-2„TJ(£>2 _ D 2) 
i<]_ , 

xh [tr (E _1 V'D 2 V + SI) - 2 tr /i'S^V'DH] (HdH') 



dF v ,D,H(V,D,H) = — ^ (VdV')(dD) 

S 2 



Expanding in power series 

2 -„| D |JV-l+*-2„TJ( D 2 _ D 2 



dF v ,D,H(V,D,H) = — ^ (VdV')(dD) 

|S| 2 

OO 

J2l h(t) [ tr (S _1 V / D a V + «)] [tr(-2tr At 'S- 1 V'DH)] t (HdH') 



x 

t=0 



From Lemma [ 

/ [tr (-2tr M 'E^V'DH)] 2t (HdH') = £ if^C. (M-WV) . 

(-) 4* 

Observing that %M , = — t , the marginal joint density of V, D is given by 

n^f | D |iv-i+K-2 n TJ (£) 2 _ D 2) 

'//v.niV.D: 



r„(f)|s| 



~ ft( 2t ) rtr(s-VD 2 v + n)l , , , a w /w N 



X 

t=0 



Now, note that D contains n coordinates for which, under this method, its corresponding 
joint density is termed cone density (or size-and-shape density) . Then we have the first main 
result of this section. 
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Theorem 2.2. The cone density is given by 



2"7T » ID^-^- 2 " Yl (Df - D 
/d(D) = - 



2^ 

3 i 

i<3 



T n [\K] r„ [|(iV-l)] |S|f 

fe(^+')(trn)A^C (D 2 ) (S^nST 1 ) 

x fc,lk *H\K) K C,(I N _ X ) ' U 

where the notation of the sum operators, CV S and A/ K are given in I Davit (1980), in par- 

6k Cl' K (1,1) 
ticular A ,' K : 



Proof. The joint density of V, D is 

lDl N - l + K -"]l(D*-D] 



X 



X 



EE 



M 2t ) [tr (S- 1 V'D 2 V + n)] C K (OS- 1 V'D 2 V) 



Assuming that ft,( 2t )(-) can be expanded in power series, 

M 2t) [trS _1 V'D 2 V + trn] - ^^^^MftrS^V'tfvf 

- £ E ^^c.( S -Wv), 

where Ce(A) is the zonal polynomial corresponding to the partition 9 = (h, ■ ■ ■ ,l a ), with 
V'/ , /, /. 

From iDavisl (|l98dh . eq. (4.13), the integration of cLFv,d(V, D) with respect to V e 
V rii Ar_i results 

/ C* K (r2S" 1 V'D 2 V) Cg (S" 1 V'D 2 V) (VdV) 

_ v ^(s-Sns-^c^ (D 2 ,D 2 ) 

And by Davisl ( 1980h . eq. (5.1) we have that 

" (W - 1+A ' ) \jy\N-X+K-2n 

i<3 



2 n n * IDI^-^" 2 ™ Yl (£>| - £ 2 ) 
/d(D) = - 



r n [\k] r„ [§(at-i)] |s|f 
fe( 2 *) (tr^A^c^p^c^Cs-Sns- 1 ) 

Now, we can derive the isotropic version of the cone density 
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Corollary 2.1. Let fcE = a 2 I, then 



/d(D) 



r„ [\K] r„ [|(iV-l)] 



1)K 



\ ^ h^jtrfl + ^trD 2 )C K (^D 2 ) g K (O) 



OC 

X 

t = K 



Proof. From Theorem 

1. n = s V' = -4 a*© V 

2. From lDavid (|l980l ). eg. (5.7), 



i 



CT^ (J 

^'"(I^n), 

2t+l \6,k 



1 \ 2t +< Aj"C^(Iiyr_i)C„(n) 



C k (Ijv-i) 

Therefore the second line of ([2]), denoted by J, it is simplified as follows: 

y> V fe< 2W) (tr»)A^C (D 2 ) (IT 1 , OS' 1 ) 



~ _ /' ri ' + ''(tvJ2i (a'';'' ') r..(D^)r,,.«. 



Note that £ 0ee . K ( A K ) Z (D 2 ) = (D 2 ) C e (D 2 ), see lDavid (Il980l) . eq. (5.10). Thus 



0,k <pe8-K 

2 



j /i( 2t +')(tr»)C K (D 2 ) C 9 (D 2 ) C K (O) 

£ tm(a 2 ) 2 ^ (±K) k C k (1n-i) ' 
- " fr( 2 ^)(trft)C K (D 2 )g K (ft) v 

» " /^+0(trft)C K (D 2 )q.(ft) , 

r&^'(^) 2t+i (K)«^(iiv-i) 11 

Now, observe that h(v) = h ^ (v — a) 1 , with a — trfi, u = trfi + ^j-trD 2 , and 

= M 2t )(v). Thus 

= A fe (2t) (tr ft + ^ tr D 2 )C K (D 2 ) C K (SI) 
~ ftpt) ( tr ft + ^ tr D 2 )q (^D 2 ) C K (O) 
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NowletbeW = D/r, r = ||D|| = ||V'DH|| = ||Y|| and noting that if D = diagpx, D n ) 
we define vecp(D) = (Dx, . . . , D n ), then 

/ D x \ I Dx/r \ 

vecp(D) = : , implies that vccp(W) = : = vecp — , 
\D n ) \ D n /r ) 

thus 

m 

(<fW(u)) = r m Y[ sin" 1 " 1 0<(du) A tfr 

i=l 

= r m J(u)(rfu) Arfr, 

with m = n — 1, u = (dx, ■ ■ ■ , 9 m ) . 

The shape density under Le and Kendall's approach it is known as disk density. 

Theorem 2.3. The disk density is given by 

n 

JJ j.W-1+Jf-an JJ (j.a _ ,.2) J(u) 
z=l i<j 

r„ [%k] |s|f 

1 /■ c K (nE- x v'w 2 v) 



/w(W) 



9 

< / s ^ N+K - n - 1 ^+ 2t - 1 h^ (s 2 + an) (d s ), (4) 

j o 



where the number of landmarks N are selected in such way that n(N+K-n-l) +t 

is a positive 

integer, then 9 = (/,, . . . ,l a ) is a partition of n ( N + K - n - 1 ) _|_^ an j Y^ =l U = n ( N+K ^~ n 1 ) _|_^ 
Proof. From Theorem 1 2. II 



dF V)D (V,D) = -gLg (dD)(VdV') 



x 



EE 



r„ [\k] |s|f 

/t( 2 *) [tr (E^V'D^)] C K (nE- x V'D 2 V) 



Let W = diag(/J, . . . , l* n ), I* = ^,r= ||D|| = ||Y||, then 



tt^ \rW\ N - 1+K - 2n r m J(u) Y[ r 2 (l* 2 - I) 
l<l 

?n[\K] |S|f 

h {2t) [ r 2 tr ( S -i V 'w 2 v)] c K (r a nir 1 v / w a v) 



If) 

dF v ,w(V,W) = f (dW)(VdV') 



EE' /'(I/O 

t=0 K 6 " V2^; 



Note that 

1 IrWI^ -1 ^' -2 " = r n{N-l+K-2n) 
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2. Ip (if If) = (r 2 )^ II (if If). 

i<j i<j 

3. c K (^ns-^'w^) = r 2t C K (nE _1 V'W 2 V). 

Collection powers of r and denning r = r , with dr = — r and 

(trS- 1 V'W 2 V) 5 (trS- 1 V'W 2 V)' 



roc 

/ r n(N+K-n-l)+2t-l h (2t) ( r 2 tr £-1 V ' W 2 V ) dr 
JO 

/ \ n(N+K-n-l)+2t-l 

= r f r) ^ (2t) (« 2 ) 



(AT+if-n-l)+2t-l 



(trS- 1 V'W 2 V) 5 / (trS" 1 V'W 2 V) 5 

(tr S" 1 V'W 2 V) 2 +J / s n(N+K-n-l)+2t-l h (2t) / 2) ds 



Thus the marginal density of dFw (W) is given by 

n^f | W |JV-l+«:-2n J (u) TJ _ ,.2) 



r„[f]|E|f S « * ! (^) K 

x / C K (nS- 1 V'W 2 V) (trE^V'W^V) 2 (VdV') 

X / s n(N+K-n-l) + 2t-l h (2t) ^2) dg 
JO 

Now, let be 

r c K (ns^vw^v) (VdV) 

^V6V„, w _i (trS-^'W^) 2 + 

Uev^-r ]T C e (5T 1 V'W 2 V) 

e 

where the number of landmarks N are selected in such way that "( JV +-^~"~ 1 ) _|_^ j s a positive 

integer, then = (/;,..., Z a ) is a partition of n(JV+ ^~"~ 1) +t, and £? =1 = n(N+K ~ n ~ 1] +t. 
Then we obtain the desired result. □ 
The isotropic case of the disk distribution follows 

Theorem 2.4. The isotropic disk density is given by 

^niH+K^ z . iV _ 1 + ^_ 2n -Q _ ^ J(u) 

/ W (W) = — 1< \n-i )K 

T n [\K]T n [£?] {a^-^~ 
C - (\ w2 ) C « ( n ) 



(5) 



v ^ ^ 

4- t!(K) K C K (I W -i) 
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Proof. The result is obtained from (Q} taking X = er 2 I and observing that |S| = |c 2 I| 

( CT 2)(iV-l). 

trS- 1 V'W 2 V = -^trV'W 2 V = ^trW 2 W = ^trW 2 = \, 
a 1 a z a 1 a 1 

recalling that ||W|| = 1. Hence 

c K (rjs^ 1 vw 2 v) (VdV) 



J 



n(Jf+K-n-l) . . 



vey„ lW -i (trE _x V'W 2 V 

2 \n(AT+Jf-n-l)+i 



/ c K ( ^nvw 2 v) (vdv') 



/ 2 sn(iV+if-n-l)+« 2n7r „ (A r_ 1)/2 C* K f^^J C« (W 2 ) 

r n [|(jv-i)] c^Iat-i) ' 

Finally, ([5]) is obtained making the change of variable s = r/cr with ds = dr/a in (U) and 
observing that C K ( \n) C K (W 2 ) = C K (n) C K (%W 2 



Alternatively, let If = A A, trD 2 = ^™ =1 £> 2 = r 2 in ©, therefore 



*2\ 
3 

i<j 



/w(W) = ,, 

r„ [\k] r„ [i(iv-i)] (a 2 )^ 
~ fo( 2t > (tr n + gj c K (g w 2 ) q (O) 

Observe that 

1 J-J" j*^./V-l+i<"-2ri _ r n(JV-l+K-2n) TJ» JV-l+K-2n 

2 - n> 2 (^* 2 - i f) = r " (n ~ i} n (^* 2 - i f) 

i<j i<j 

3. C K (^W 2 ) =r 2t C K (^W 2 ). 
Collecting powers of r we have 

J™ r n(N+K-n-l)+2t-l h (2t) ^ O + ^ dr. □ 

3 Central Case 

We obtain the central cases of the cone and the disk densities. 
Corollary 3.1. The central cone density is given by 

2 n \n\ N - 1+K - 2n T\ (Df - D?) 

f (m _ tj f r h"'(0)C,(E- 1 )C,(D») 

/D( ] «- sS! ^T n [iK]r n [XN-i)]m*hr flG «c^-o ' {) 
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Proof. Start with 

dF Y (Y) = [trE^YY'] (dY). 

The joint density of V, D, H is 

J] {D\ - D>) (rfD) 

dFV^(V,D,H) = ^Z^Z^^ [trSTWv] (HdH')(VdV). 
Recalling that 



/ (H'dH) = 

./HeV„,K 



r„ [ijr] 



we have that 

tt^ | D |*r-i+if-2n H (pi _ (dD) 

dF v , D (V,D) = n-^-x h [trE^V'D'V] (VdV) 

1 „ L 5 AJ |E| 2 

Then the required density is given by 



| D |*-i+*-a»IJ(D?-.D») 
= _^ / ft [trSTWv] (VdV). 

Integrating with respect to V 

/ h [tr IT 1 V'D 2 V] (VdV') 

JV n ,N-l 

= f jr (trE^V'D^)' (VdV) 

•/V„,jv-1 ; =0 

= EE^/ C * (ST 1 V'D 2 V) (VdV) 

;=0 ' "'Vn.JV-l 



~ ^(0(0) 2"7r^C (ST 1 ) (D 2 ) 

~ v ft«)(o) c e (sr 1 ) c e (d 2 ) 



Then, 



/d(D) 



i<3_ 

r n [\k] r„ [Un-1)] |s|f 



" V / 1 W(Q)C 9 (S- 1 )C 9 (D 2 ) 



Alternatively, from Theorem 12.21 if we take 
1. /i< 2 *+') (trJl) = ftW(O), 
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A e,K _ A e,e _ Ce(I) _ 1 
l - ^ ~^ ~ C e (I) ~ ' 

3. C (D 2 ) = G,(D 2 ), C*(I) = C e (T), 

4. (E- 1 , OS- 1 ) = (S- 1 , 0) = C e (IT 1 ) , 

the required result follows. □ 
And finally, note that the central disk density is invariant under the elliptical distribu- 
tions. 

Corollary 3.2. The central disk density is given by 

n 

ri±(n(N +K-n- 1))] n n N - i+K - 2n n of - ^ 2 ) J W 



/w(W) = 



2^^r„ [|K]|E|f 

(VdV) 



vev„,„_i (E _1 V'W 2 V) 



(7) 



where the number of landmarks N are selected in such way that n ( N+K " — JO j s a p 0S m ve 
integer, then 9 = (li, . . . ,l a ) is a partition of n ( N+K anc i \ i = "(N+if-n-i) 

a 

where 9 = (h, . . . ,l a ) is a partition of the positive integer n(N + K — n — l)/2, / ^h = 

i=l 

n{N + K-n-l)/2. 

Proof. From Theorem^ taking t = 0, il = 0, h^°\-) = h(-) and recalling that 

7T 2 

the result is obtained. □ 

4 Example: Mouse Vertebra 

As the reader can check the general cone and disk densities are given in terms of invariant 
polynomials, so at this time no inference can be performed, except if the series are truncated 
in the first few terms. However, there is a way to work with an exact density, it is, when we 
assume an isotropic model. 

Consider the isotropic elliptical disk density of theorem 12.41 



/w(W) 



2 n n ^^ Y[ I* N ~ 1+K - 2n J] (If - If) J(U) 



r n [\k] r„ [|(iv-i)] {a*) 2 ^- 



x 



c K (^w 2 ) c K (n) 



^t\(lK) K C K (lN-i) 

r n(N+K-n-l)+2t-l h {2t) ^ n + ^ ^ 
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and the generator for the subfamily Kotz 



Ky) 



R 1 



..T-l . 



-Ry 



77 K(N-1)/2 T (T - 1 



A'(jV-l) 
2 



with derivative 



dy 



-_y T - x e~ R y = {-R) k y T - x e 



R y(i 



E 



i=0 



{-Ry)' 



see ICaro-Lopera et al ( 2009t ) for other families (Pearson VII, Bessel, general Kotz, Jensen 
Logistic) and their derivatives. 

Now, we contrast three models next, the classical Gaussian (T = 1, R = |) and two 
non normal (T = 2, R = | and T = 3, R — ^) via the modified BIC criterion, they will be 
applied to the data of two group s (sma ll and large) of mouse vertebra, and experiment very 
detailed in lDrvden and Mardial (|l998l) . 

The isotropic Gaussian disk density is given by 



2^{-2-M+n+Kn-n 2 +nN) J~J (* N-l + K-2n J~J ^*2 ^*2~j J^ u ~j 



/w(W) 



7r ^M+n-Kn-nN) a M-n(-l+K-n+N) Tn [ljf] [1 (jy _ ]_) 



x ctr 



~ r [1 + K - n + TV)) + t] v ^ (W 2 ) C K (f^ 



where M = K(N - 1) and n = min{(7V - 1), K} . 

The Kotz disk density when T — 2 and R—\ follows after some tedious simplification, 
and it is given by 



/w(W) = — 



2^(~M+n+Kn-7i 2 +nN) J~J [* N- l + K-2n J~J ^*2 _ ^*2~j 



7r i(M+n-Kn-nJV) (7 M-n(-l+K-n+JV) rr[ ^ [1 (jy _ !)] M 



x etr — 



rr2 J +! 



2a 2 



(=0 



tr 



2cr 2 



2t r 



: (n(-l + X - n + N)) + t + 1 



'-(n(-l + K -n + N))+t 



c K (w 2 ) C K (0) 
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Finally, the corresponding density for the Kotz model T = 3, is obtained as: 



/w(W) = 
„ J ( u ) 



2^(-2-M+n+Kn-n 2 +nN) J~J £* N-l+K-2n J~J ^*2 _ ^*2~j 

i=l 

n i(M+ n -Kn-nN) a M-n(-l+K-n+N) Tn [ljf] ^ [l(JV-l)] 



M{M 

OO 



i! 

t=o 

+4 I tr 



2, 
-8t 

M> 

2a 2 



■etr 



m 1 



4± r 



16t tr 
1 



MM 

2a 2 



4tr 



^2a 2 



(n(-l + X - n + iV)) + t + 1 



-4r 



(n(-l + K-n + N)) +t + 2 



E 



c K fw 2 ) a. 



MM 

2cr 2 



The likelihood based on exact densities require the computation of the above series, a 
carefully comparison with the known hypergeometric of two matrix argument indicates 
that these distribut i ons ca n be obtained by a suitable modification of the algorithms of 
iKoev and Edelman I (|2006h . 

In order to decide which the elliptical model is the best one, different criteria have been 
employed fo r the model selection. We shall consider a modification o f the BIC statis tic as 
discussed in lYang and Yane ( 2007t) . and which was first achieved by Rissanen (1978) in a 
coding theory framework. The modified BIC is given by: 

BIC* = -2£(£, a 2 , h) + n p (\og(n + 2) - log 24), 

where £(/i, cx 2 , h) is the maximum of the log- likelihood function, n is the sample size and n p 
is the number of pa rameters to be estimated for ea ch particular sh ape density. 

As proposed by iKass and Rafterv" ( 1995t ) and Rafterv ( 1995 ) , the following selection 
criteria have been employed for the model selection. 



Table 1: Grades of evidence corresponding to values of the BIC* difference. 



BIC* difference 


Evidence 


0-2 


Weak 


2-6 


Positive 


6-10 


Strong 


> 10 


Very strong 



The maximum likelihood estimators for location parameters associated with the small 
and large groups arc summarised in the following table: 

According to the modified BIC criterion, the Kotz model with parameters T = 2, R = h 
and s = 1 is the most appropriate, among the three elliptical densities selected, for modeling 
the data. There is a very strong difference between the non normal and the classical Gaussian 
model in this experiment. 

Let n l and fi 2 be the mean disk of the small and large groups, respectively. We test 
equal mean shape under the best model, and the likelihood ratio (based on —2 log A « Xio) 
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Table 2: Maximum likelihood estimators. 



Group 


BIG* 


fill 


Ml2 


fJ-21 


A*22 


fJ-31 


M32 




G 

K : T = 2 
















K:T=3 
















-518.3844 


2.3931 


-27.7470 


34.2670 


-16.0397 


-9.9515 


0.4755 


Small 


-536.1662 


-16.3534 


-75.7142 


40.2639 


-3.9264 


39.6012 


-15.9524 




-116.7475 


-1.1862 


-19.3374 


41.8450 


-15.4168 


44.9495 


0.3280 




-519.1378 


-24.8027 


-63.4131 


35.8304 


-9.0617 


31.9550 


-15.0595 


Large 


-536.9652 


-19.3968 


-69.2566 


41.9566 


-1.7093 


32.1251 


-20.2930 




-115.5510 


-27.1256 


-61.2295 


36.5588 


-10.8047 


73.2563 


-6.8034 













fJ-41 


fJ-42 




M52 


a 2 


17.2967 


24.2432 


14.8582 


13.5573 


33.9149 


16.0621 


15.5973 


-52.0684 


13.4224 


106.7591 


21.6988 


18.9804 


-5.8144 


10.1334 


26.7961 


11.4842 


14.7148 


-44.1020 


21.2210 


94.1200 


13.2026 


13.8409 


-48.3278 


16.6046 


94.2873 


16.8063 


17.3769 


-46.5415 


23.0966 


71.1297 



for the test H : fi 1 = fi 2 vs Ha '■ Mi 7^ A*2) provides the p- value 0.92, which means that 
there extremely evidence that the mean shapes of the two groups are equal. 

It is important to note that the general densities derived here apply to any elliptical 
model; some classical elliptical densities as Kotz, Pearson II and VII, Bessel, Jensen-logistic, 
can be obtained explicitly and applied, however they demand the computation of the the 
A:-th derivative of the generator elliptical function h(-). This is n ot a trivial fact, but for th e 
above mentioned families, the required formulae are available in Caro-Lopera et al ( 2009h . 



However, isotropic densities are more tractable because they are expanded in terms of zonal 
polynomials, instead of non isotropic distributions which require some additional conditions 
on the number of landmarks in order to obtain a know integral expanded in terms of invariant 
polynomials. The general densities expanded in terms of invariant polynomials, seem non 
computable at this date for large degrees. 
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